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20.  Abstract  (Continued) 

from  which  the  dynamic  characteristics  of  a heated  counter  flowing  jet  in  the 
absence  of  buoyancy  can  be  calculated.  The  present  investigation  is  concerned 
with  the  effect  of  buoyancy  upon  the  motion  of  a counterflowing  jet. 

The  lower  portion  of  the  trajectory,  which  has  been  calculated  by  the 
present  model,  is  in  fair  to  good  agreement  with  the  corresponding  experi- 
mental curve,  the  calculated  curve  tending  to  be  somewhat  higher  than  that 
obtained  experimentally.  The  calculated  upper  part  of  the  trajectory,  obtained 
from  a model  which  gives  the  deflection  of  a jet  in  a crosswind,  is  in  good 
agreement  with  experiment.^^ 

The  present  model  yields  a scaling  law  for  a counterflowing  jet  which  indi- 
cates that  the  scaling  depends  principally  upon  the  Froude  number  defined  by 
the  initial  jet  velocity  and  excess  temperature  and  very  weakly  upon  the  initial 
jet  temperature.  This  result  is  essentially  the  same  as  that  obtained  for  the 
coflowing  jet  for  the  case  of  small  values  of  windspeed  relative  to  initial  jet 
velocity. 
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Calculation  of  the  Buoyant 
Motion  of  a Turbulent  Planar 
Heated  Jet  in  an  Opposing  Air  Stream 


1.  INTRODUCTION 

As  part  of  the  development  of  an  operational  Warm  Fog  Dispersal  System 
(WFDS),  experimental  and  theoretical  studies  have  been  made  of  the  characteris- 
tics of  ground  based  heated  jets-  for  various  combinations  of  heat  and  thrust  under 
different  wind  conditions.  The  dynamic  characteristics  of  a nonbuoyant  coflowing 
jet,  that  is,  jet  in  the  same  direction  as  the  wind,  are  well  known  and  presented  in 
detail  by  Abramovich. 1 The  method  of  calculating  the  buoyant  motion  of  a heated 

O 

submerged  jet,  that  is,  no  wind,  (Abramovich  ) can  be  extended  in  a straightforward 
manner  to  the  case  of  the  coflowing  jet,  planar  or  round  (Klein  and  Kunkel  ’ ).  For 
simplicity,  the  attachment  of  the  jet  to  the  ground  (ground  effect)  was  neglected  in 


(Received  for  publication  22  March  1978) 

1.  Abramovich,  G.  N.  (1963)  The  Theory  of  Turbulent  Jets,  The  MIT  Press, 

Cambridge,  Mass.,  Chapters  4 and  5. 

2.  Abramovich,  G.  N.  (1963)  The  Theory  of  Turbulent  Jets.  The  MIT  Press, 

Cambridge,  Mass.,  Chapter  1 i,  pp  580-585. 

3.  Klein,  M.M.,  and  Kunkel,  B.A.  (1975)  Interaction  of  a Buoyant  Turbulent 

Planar  Jet  With  a Coflowing  Wind.  AFCRL.-TR-75-0368. 


4 Klein  M.M.,  and  Kunkel,  B.A.  (1975)  Interaction  of  a Buoyant  Turbulent 
Round  Jet  With  a Coflowtng  Wind,  AFCRL-TR-75-6581" 


these  investigations.  A method  of  taking  into  account  the  ground  effect  for  a co- 

5 

flowing  jet  has  been  developed  by  Klein. 

The  situation  for  the  nonbuoyant  counter  flowing  jet,  that  is,  jet  and  wind 

g 

directions  opposite,  as  presented  by  Abramovich,  is  considerably  less  satisfac- 
tory both  with  regard  to  theory  and  experiment.  Here  the  calculations,  which  must 
account  for  the  regions  of  counterflow,  are  quite  specialized  and  not  easily  adapted 

to  the  buoyant  jet.  Recently,  however,  a simplified  model  for  a counterflowing 

7 

round  jet  has  been  developed  by  Sekundov,  which  may  be  extended  in  a straight- 
forward manner  to  take  account  of  buoyancy  and  obtain  the  jet  trajectories.  An 
important  feature  of  the  Sekundov  model  is  the  use  of  a finite  wall  to  help  simplify 
thq  equations  of  motion.  The  results  for  an  open  jet  are  then  obtained  by  making 
the  wall  arbitrarily  large.  The  Sekundov  model,  which  is  limited  to  the  case  of 
an  unheated  incompressible  jet,  has  been  extended  to  take  account  of  heat  addition 

g 

and  density  variation,  while  neglecting  the  buoyant  motion  (Klein  ).  The  present 
investigation  is  concerned  with  the  buoyant  motion  of  the  counter  flowing  jet.  Since 
the  jets  in  the  warm  fog  dispersal  system  merge  a short  distance  downstream  of 

Q 

the  jet  nozzles,  the  investigation,  parallel  to  that  reference  8,  has  been  confined 
to  the  planar  case. 

g 

As  in  the  case  of  the  coflowing  jet,  experiments  show  that  the  counterflowing 
jet  also  remains  attached  to  the  ground  for  some  distance  downstream  of  the  nozzle, 
resulting  in  a delayed  lift-off  point.  A method  of  analysis  similar  to  that  of  refer- 
ence 5 was  utilized  to  determine  the  point  of  lift-off. 


2.  JET  GEOMETRY 


A schematic  sketch  of  the  flow  pattern  for  a planar  jet  is  shown  in  Figure  la 
while  the  corresponding  trajectory  of  the  jet  centerline  is  given  in  Figure  lb. 
Here  u and  T denote  velocity  and  temperature,  while  the  subscripts  o,  m and  a 
designate  initial,  axial  and  ambient  values.  In  the  initial  section,  where  the 


Klein,  M.M.  (1977)  A Method  for  Determining  the  Point  of  Lift-Off  and 
Modified  Trajectory  of  a Ground-Based  Heated  Turbulent  Planar  Jet~ 
in  a Coflowing  Wind,  AFGL-TR-77-d()33. 

Abramovich,  G.N.  (1963)  The  Theory  of  Turbulent  Jets,  The  MIT  Press, 
Cambridge,  Mass. , Chapters  4 and  9. 

Sekundov,  A.N.  (1969)  The  Propagation  of  a Turbulent  Jet  in  an  Opposing 
Stream,  in  Turbulent  Jets  of  Air,  Plasma  and  Real  Gas,  Consultants 
Bureau,  New  York. 

Klein,  M.M.  (1977)  Interaction  of  a Turbulent  Planar  Heated  Jet  With  a 
Counterflowing  Wind,  AFGL-TR-77-0214. 

Kunkel,  B.A.  (1975)  Heat  and  Thrust  Requirements  of  a Thermal  Fog  Dist 
System,  AFCRL-TR-75-0472. 


boundary  layer  has  not  yet  reached  the  center  of  the  jet,  the  surface  y^x)  separates 
the  regions  of  perturbed  and  unperturbed  flow.  In  the  main  section  the  boundary 
layer  has  reached  the  jet  axis  and  the  velocity  profile  remains  constant  in  form. 

The  surface  y2(x)  indicates  where  the  jet  velocity  has  reversed  its  direction,  while 
y3(x)  designates  the  streamline  which  separates  the  regimes  of  perturbed  and 
unperturbed  flow.  Because  of  the  constricting  effect  of  the  walls  the  flow  conditions 
in  the  region  above  y3>  designated  by  the  subscript  c,  are  almost  but  not  quite 
ambient. 

Experimental  work  indicates  that  the  jet  flow  can  be  separated  into  two  regimes, 
Zone  I and  Zone  II.  In  Zone  I,  consisting  of  the  initial  section  and  the  main  section, 
the  flow  pattern  is  qualitatively  similar  to  that  obtained  for  a submerged  jet,  that 
is,  the  pressure  is  almost  constant  and  the  surface  y2  grows  at  a constant  rate. 

Here  the  flow  pattern  can  be  accurately  calculated  within  the  framework  of  the 
model. 

In  Zone  II  the  pressure  rises  very  rapidly  and  the  surface  y^  decreases  to  zero 
as  the  final  retardation  of  the  flow  takes  place.  Although  an  accurate  evaluation  of 
the  flow  conditions  is  not  possible  for  this  regime,  interpolative  procedures  may 
be  used  since  the  extent  of  Zone  II  can  be  calculated  and  conditions  at  its  end  are 
known.' 


3.  ANALYSIS  FOR  ZONE  I 

3.1  Axial  Distributions  for  Counterflowing  Jet 

The  axial  velocity  um  and  the  axial  temperature  excess  ATm  in  Zone  I are 
given  as  functions  of  jet  distance  x by  (see  reference  8) 


10 


r 


where  AT  is  the  temperature  excess  over  ambient,  m is  the  dimensionless  value 
m 

of  the  initial  jet  velocity  uq  referred  to  ambient,  and  the  subscripts  H and  1 refer 
to  the  ends  of  the  initial  and  main  sections.  As  in  reference  8,  velocity,  density, 
temperature  and  length  are  taken  as  dimensionless,  referring  to  initial  jet  half 
width  yQ  for  length  and  to  ambient  values  for  the  other  quantities. 

The  analysis  presented  in  reference  8 shows  that  the  temperature  variation 
given  by  Eq.  (2)  is  somewhat  slower  in  the  outer  region  of  the  jet  than  that  given 
by  experiment.  We  shall,  therefore,  utilize  in  place  of  Eq.  (2)  the  linear  form 


AT 


m 


AT 


1 + 

m 

1 + m 


(3) 


which  is  close  to  Eq.  (2)  in  the  initial  region  of  the  jet  but  gives  a slightly  greater 

decrease  in  AT„  in  the  outer  region, 
m 

The  surface  of  zero  velocity,  y2(x),  and  the  outer  boundary  of  the  jet,  y3(x), 
are  obtained  from 


y2 


(4) 


N - 1 =ri-pm  « (2u  + 3) 

a m m m 
o 


(51 


1 L+PJL\  1 f , , 1 + ATm 
"1  = 2 V Pm)  '2  1 + AT2 


(6) 


ao  = i(1  + p2)  = 2 ( 1 + 1 + AT, 


(7) 


AT 

AT  = — 

1 + u, 


m 


(8) 


y o 

where  N = — — , c is  a growth  or  mixing  coefficient  having  the  value  0.  22  for  the  main 
y2 

region,  p is  density,  and  the  subscript  2 indicates  values  at  the  surface  y2-  Be- 
cause of  the  delayed  lift-off  point  and  the  initial  rapid  decrease  in  temperature,  the 
quantities  ATm  and  AT2  are  small  compared  to  unity.  It  will,  therefore,  be  con- 
venient to  develop  aj,  aQ  and  pm  to  2nd  order  in  ATm  to  yield 


al  = 1 + T71  *Tm 


. (iLll)  at2 
0 z m 
2v 


(9) 
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(10) 


1 12 
a = 1 + tt-  AT  - -V  AT*L 
o 2v  m 2V^  m 


p = 1 - AT  + AT 
m m 


N - 1 = (2v  + l)(v  - 1) 


AT, 


1 -- 


M) 


9 2/  AT2 
2 v m 


(ID 


(12) 


where  v = 1 + u 


m" 


3.2  Calculation  of  the  Buoyant  Force 

The  buoyant  force  per  unit  length,  B,  is  obtained  from 


B ■ 2gps 


P 

- p)  dy  + / (1 

^2 


- p)  dy 


(13) 


where  y is  the  \ ertical  coordinate.  In  view  of  the  small  variation  of  pressure  in 
the  main  region  we  may  write 


J_  1 

p T 1 + AT 


(14) 


and  express  Eq.  (13)  as 


B 


2gpa 


AT 

1 + AT 


dy  + 


AT 

1 + AT 


(15) 


Since  AT  is  small  compared  to  unity  at  and  beyond  the  lift-off  point,  we  may 
develop  the  integrand  in  Eq.  (15)  as  a power  series  in  AT  to  yield 


B 2gpayo(I1+I2) 

12 


I 


x = f AT(1  - AT  + AT2  - AT3  + . . . ) dy 


(16) 


(17) 
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AT(1  - AT  + AT2  AT3  + . . . ) dy 


(18) 


To  evaluate  1^  and  I2,  and  for  subsequent  analysis,  we  require  the  velocity  and 
temperature  profiles.  Utilizing  the  linear  profiles  of  reference  8,  for  the  region 

0 < y ^ y2: 

= i -X_  (19) 


AT  u + 1 , y 

ATm  " um+l  " " **2 


a = 1 + ■ 


AT  1 + u 
m m 


For  the  region  y2  ? y i y3: 


y - y2 
y3 


1 - u = 1 - 


y - y2 
y3  ^2 


we  obtain  for  Ij  and  I2, 


'l  ' *2  iTm  [>  -Ti  ' iTm  (‘  ' T 4 J?)  4 (‘  - h*  ? - ^l)]  1241 

AT,  / , AT?  \ 

I2  = y2<N  ' » ~2  V,1  ‘ f ^T2  + ~2  ) * (25) 
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Employing  Eqs.  (8)  and  (12)  for  AT,,  and  N - 1,  and  expressing  a in  terms  of  v by 


1 -I 

v 


(26) 


allows  us  to  write  Eqs.  (24)  and  (25)  in  the  form 


*1  = y2  ATn 


h = y2  ATr 


HI+U  -T8  (l  + 7+^)+^ 


V V 


(27) 


(*♦*) 


+ ATm  (1+TT7  + 72 


(28) 


Addition  of  Eqs.  (27)  and  (28)  then  yields  for  the  buoyant  force  in  Zone  I. 

B = 2gPayoy2[^Tm-  (?+|'41f)  ATL+  (|  + T2  + £h-^2) 

*Tm]  • 


(29) 


3.3  Trajectory  Analysis 

The  vertical  velocity  is  determined  from  the  equation 

^(^V)  = B (30) 

where  <t>  is  the  mass  flux,  through  the  cross  section  of  the  boundary  layer  and  given 

by 


<t>  = <t>1  + <t>2  (3D 

with  i/ij  the  mass  flux  through  that  portion  of  the  boundary  layer  below  y^  and  <t> 2 
the  mass  flux  through  the  part  above  yg.  We  note  that,  although  4>j  and  <t>2  are  in 
opposite  directions,  <l>  is  given  by  the  numerical  sum  of  ^ and  <t>2-  This  is  due  to 
the  fact  that  the  equation  for  the  vertical  velocity  is  unaffected  by  a change  in 
direction  of  the  flux,  that  is,  the  vertical  velocity  increases  positively  to  the  right 
in  either  case. 
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The  quantities  ^ and  are  obtained  from 


r 

),  = 2 p uy  I pu 
1 Ka  a Jo  J K 


dy 


4 = Pa  ua  al  pm  um  y2  yo 


(32) 


<t>2  = 2Pj 


Uayo|  P U 


dy 


^2  = Pa  ua  ao  <N  - 1}  y2  yo  = pa  ua  al  Pm  um(2um  + 3)  y2  yo 


(33) 


= 2p  u a,  p u (2  + u ) y.  y 
^a  a 1 Mm  m'  m J2  Jo 


(34) 


in  which  we  have  used  Eq.  (5)  for  N - 1. 

The  slope  of  the  center  line  of  the  jet  is  given  by 


dy  . <f>V 


M 


(35) 


where  P is  the  momentum  flux  through  the  boundary  layer  and,  analogous  to  the 
mass  flux  <j>,  given  by 


P = Pi  + P2 


(36) 


*2 


Pl  = 2PaUayo  / PU 


dy 


P1  = |paUaalPmUm  y2  yo 


(37) 


'of  P 


P2 = 2pa  ya  / pu2  dy 


P2  = I Pa  ua  ao(N  - y2  y0  = I Pa  ua  al  Pm  um(2um  + 3)  y2  yo 


(38) 
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(39) 


P = 2?a  ua  al  Pm  Vum  + 1}  • 

Since  the  point  of  lift-off  x is  well  downstream  of  the  nozzle,  the  axial  veloc- 

a ‘ 

ity  and  temperature  have  dropped  considerably  below  their  initial  values  and  vary 
slowly  beyond  x . An  inspection  of  the  results  of  reference  8 shows  that  the  veloc- 
ity um  decreases  almost  linearly  with  distance  x beyond  the  lift-off  point.  There- 
fore, for  purposes  of  integration,  it  is  convenient  to  replace  the  cumbersome 
equation  (1)  by  the  simple  linear  from 


u 


m 


u - (u 
a a 


(40) 


Utilizing  Eqs.  (3)  and  (4)  for  AT  and  y„  we  may  write  the  integrated  form  of  the 

n m 6 

buoyancy  equation  to  order  AT^  in  the  form 

9 AT  f AT  1 

♦V  = *aVa  + 2gpay2CrJE  (41) 

where 


(42) 


(43) 


Numerical  checks  show  that  the  J integral  contributes  very  little  to  the  vertical  ve- 
locity V and  will  therefore  be  neglected.  The  I integral  is  easily  evaluated  to  yield 


16 


where 


3 = u + 1 + (u  - 1)  t 
a a a 


t = 


x.  - x 
1 a 


(45) 

(46) 


Making  use  of  Eq.  (34)  for  <t>,  we  may  write  the  vertical  velocity  equation  in 
the  form 


x (v  - 1)  f V , _ 

V _ a a q q *■  m 


uo  ' x V2  - 1 


__  ...  1 J_(I  J . 

f u K 1 + m 2 . x o’ 

O V - I 


(47) 


where 


K = 


y0gATo 


(48) 


is  the  Froude  number  for  the  flow,  and  f is  the  value  of  alPrn  in  Eq.  (34),  and 
given  to  first  order  in  AT^/(1  + m)  by  (see  Eqs.  (9)  and  (11)) 


f = 1 - 


1 + v *To 
2 1 + m 


(49) 


If  we  employ  Eq.  (41)  for  the  vertical  velocity,  the  integrated  form  of  the  tra- 
jectory Eq.  (35)  may  be  written  as 


(50) 


Using  Eq.  (39)  for  P,  the  integrals  in  Eq.  (50)  are  easily  evaluated  to  yield 


««»  - K.>  4 1^5  «k2  ‘ Ki> -k  ■„<*»  ' Ki>  l51> 


17 


where 


ko-^t  / 


dt 


t(v  - l)  ' 3 - l 


^-T  in  (tH) 


(52) 


Ki  = f/tV  =Ttn  (v) 


(53) 


k2  j 


XdL 


t(v  - 1) 


dt 


K„  = 


(3  - l) 

(v  - 2)2 

O' 


. t 1 \ a2  . , 53  + 3 .2  (va  " 2) 

a2<n  (v-l)'v  - 2 1 24  1 " 12  1 


(54) 


K„ 


t 

f L^t 

J tv 


*3  “ 


3s 


12(v  - 2) 

a 


in  i 


CO  9 (v  - 2) 

_ _ * + JLiLt*  2 

v 12(v  - 2)  12  12 


r 


(55) 


= 11  _ (53  + 3X3  - 1) 


2 " 2 


12 


(56) 


and  we  have  neglected  the  small  term  in  ATq/(1  + m)  contributed  by  f. 


3.4  Dependence  of  Vertical  Velocity  and  Trajectory  Upon  Parameters 

The  foregoing  results  give  the  general  dependence  of  the  vertical  velocity  and 
trajectory  upon  the  jet  parameters  and  position.  Thus,  from  Eq.  (47)  the  vertical 
velocity  is  proportional  to  1/K  while  Eq.  (51)  indicates  the  trajectory  y is  pro- 
portional to  m/K  . This  explicit  dependence  upon  m is  due  to  the  separation  of  m 

3 

and  x in  the  velocity  Eq.  (1).  It  is  not  present  in  the  trajectory  for  a coflowing  jet 
where  the  dependence  of  vertical  velocity  upon  windspeed  is  far  more  complex. 

An  examination  of  Eqs.  (47)  and  (51)  shows  that  the  dependence  of  V/uq  and  y 
upon  position  x is  more  complex  than  a simple  power  law.  However,  an  analysis 
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of  numerical  solutions  of  these  equations  indicates  that,  roughly,  the  vertical 

2 

velocity  increases  less  rapidly  than  x while  the  trajectory  rises  more  rapidly 
2 

than  x . These  appears  to  be  reasonable  results  when  compared  to  the  coflowing 
_ case  where  V/uQ  ~ x>  y~x^^- 

3.5  Scaling  Law  for  Counterflowing  jet 

The  foregoing  results  may  be  utilized  to  derive  the  scaling  law  for  a counter- 
flowing jet.  An  examination  of  Eqs.  (1)  and  (2)  show  that,  for  constant  m,  the 
velocity  and  temperature  distributions  are  independent  of  length  scale.  The  right 
hand  side  of  the  trajectory  Eq.  (35)  will,  therefore,  depend  only  on  the  Froude 
number  K during  a change  of  scale  at  constant  m (note  that  the  velocity  dependence 
v ~ Tq  in  both  numerator  and  denominator  leads  to  a weak  dependence  upon  Tq 
which  may  be  neglected  here).  This  scaling  law  is  in  conformity  with  that  obtained 

3 

for  the  coflowing  jet  where  it  is  shown  that  the  trajectory  depends  principally  upon 
the  Froude  number,  the  effect  of  initial  temperature  Tq  being  negligible. 

4.  ANALYSIS  FOR  ZONE  II 
4.1  Axial  Diatributiona 

As  indicated  previously,  the  axial  velocity  and  temperature  in  Zone  II  do  not 
have  a large  variation  and  may  be  obtained  by  interpolation.  Since  um  has  the 
value  unity  at  x^  and  zero  at  the  end  of  Zone  II  we  write 


u = 1 — S (57) 

m 

where 


and  x2  designates  the  end  of  Zone  II.  The  experimental  data  indicates  that  Eq.  (3) 
for  ATm  is  fairly  accurate  in  Zone  II  and  will,  therefore,  be  retained  in  this  region. 

The  surface  y2(x)  decreases  to  zero  with  a vertical  tangent  at  x^;  accordingly, 
we  shall  utilize  the  simple  parabolic  form 

y2  = y21d  - S)1/2  (59) 

where  y2^  designates  the  value  of  y2  at  x^. 
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The  experimental  data  of  Vulis10  indicate  that  y3  is  roughly  constant  until  y2 

and  u have  decreased  to  zero.  A short  distance  further  out  ambient  conditions 
m 

have  been  attained.  For  simplicity  we  shall,  therefore,  assume  y3  is  constant  and 
take  the  end  of  Zone  II,  Xg,  as  the  termination  of  the  jet,  that  is,  the  location  at 
which  the  jet  reverses  its  direction.  Accordingly,  we  may  now  write 


*3  = *31  = y21Nl 

where  y31,  and  are  the  values  of  y3  and  N at  x^. 
is  given  by 


(60) 

From  Eqs.  (3)  and  (12), 


Nj  - 1 = 5 


1 


*To 
1 + m 


+ 2. 

2 


AT" 


(1  + m) 


(61) 


4.2  Determination  of  the  Buoyant  Force 

The  buoyancy  at  the  position  x^  is  given  by  the  integrals  Ij  and  1^  in  Eqs.  (24) 
and  (25),  but  the  values  of  y2,  y3,  v and  N are  now  given  by  the  equations  of  Section 
4.  1.  Accordingly  we  first  write  the  buoyancy  B as 


B = 2 gp. 


yl  (B 


1 + B2  + B3) 


(62) 


where 


i at 
1 o 


B1  = 2 TTIH  (y2  v + y3> 


AT 

B„  = - I 2_  [y2(v2  + v)  + y3] 


2 3(1+m)2-2 


1 


B„  = i- % [y2(v3  + v2  + v)  + y3)] 

3 4 (1  +m)3  2 J 


(63) 


(64) 


(65) 


and  then  use  (59),  (60)  and  (61)  to  cast  Bv  B 2 and  B3  into  the  more  explicit  form 

i AT  a /« 

Bi = \ its  y2i  «2  - s)(1  • s) + 61  (66) 


10.  Abramovich,  G.  N.  (1963)  The  Theory  of  Turbulent  Jets.  The  MIT  Press, 
Cambridge,  Mass.,  Chapter  1,  pp  32-3(5. 
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T 


B„ 


B 


i AT^  y21  j «2  - S)2  + (2  - S)J(1  - S)1/2  + -y  | 
3 (1  + mV  41  I » 


- -AT°-t,  y„,  | [(2  - S)3  + (2  - S)2  + (2  - S)](l  - S)1/2  +^[ 

4 (1  + mr  21  I > 


3 4 (1  + m) 

4.3  Calculation  of  the  Trajectory 

The  flux  if)  in  Zone  II  is,  from  Eqs.  (32)  and  (33)  now  given  by 
* +*2  +^3 


y?  y3 
*1  = 2^(v  - 2)  +_r 


(67) 


(68) 


-AT 


u2  = 4(1  +°m)  Iy2(v  ' 2)  + y3] 


AT. 


^3 


4(1  + m) 


° — 2 {y2I<v  ' 1Hv  + 1)  - 1]  } 


(69) 


(70) 


(71) 


(72) 


while  the  momentum  P is,  with  Eqs.  (37)  and  (38),  now  defined  by 


P = + \2  + x3 


1 9 ^3 

X1  =7y2l(v-  l)  - 11  + 2 


-AT 


X2  = 60~T^)  <y2  I(v  - 11  <v  + M + y3 } 


(73) 


(74) 


(75) 


AT 


x3 = "~°>2  {y2l(v  - 1)2(y2  + i»  - 1]  + y3> 

6(  1 + m) 


(76) 
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Using  Section  4.  1,  we  can  develop  the  flux  and  momentum  equations  to  the  form 


Pj  = -5^  [6  - S(1  - S)1/2] 

n2  - -t1  TtL  (K1  ' SH3  - s>  - 1 1 < 1 " S)1/2  + 16} 

AT2 

= Z|i  o {[(1  - S)(2  - S)2  - S]  (1  - S)1/2  + 26} 

3 2 (1  + mr 

A1  =~f1  f[(1  ‘ S)2  ' 1,(1  " S>1/2  + 6> 

X2  = ^ -^2.  (((i  . S)2(3  - S)  - 1 ] ( 1 - S)1/2  + 16} 

2 

x .111  ATo  r[(1  . s)2(2  . g)2  + (1  . s)2  - 1]  <1  - s)1/2  + 26}  . 
3 6 (1  + mr 


(77) 

(78) 

(79) 

(80) 

(81) 

(82) 


Because  of  the  dominant  effect  of  the  y3  term,  the  flux  and  momentum  terms  vary 
very  little  in  the  Zone  II  region.  We  shall,  therefore,  utilize  average  values  of 
these  quantities  in  calculating  the  vertical  velocity  and  trajectory.  Integration  of 
the  buoyancy  equation  now  yields  for  the  vertical  velocity 


V 

u 


o 


xl)  K 


1 + m 


2 L 


(83) 


where 


L = 


1 - r3/2 
3 


+ 


+ 3(1  - r) 


u = 6 - S(1  - S)1/2 
r = 1 - S 

in  which  we  have  neglected  the  small  contribution  of  the  ATq  terms. 


(84) 

(85) 

(86) 
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Integration  of  the  trajectory  equation  yields 


3 „ V1  6 


y = yl  m—  X (x’xl)  + r 


m 


+ m 


M 


= M s . 2_(  S/2j 
15s  15  ' ' 35  u 


A = 6 + [(1  - S)  - 1]  (1  - S) 


1/2 


ri(x2-V2m 

r7/2).3  2 

2 


(87) 

(88) 

(89) 


5.  CALCULATION  OF  UPPER  PORTION  OF  TRAJECTORY 


The  calculated  results  show  that  the  slope  of  the  trajectory  at  the  end  of  Zone  I, 
x, , is  generally  near  unity  and  increasing  rapidly  as  it  moves  toward  a vertical 
position  at  Xg.  During  this  interval  the  character  of  the  motion  has  altered  from 
one  in  which  the  change  of  direction  is  due  principally  to  buoyancy,  to  one  in  which 
it  stems  primarily  from  deflection  due  to  wind  pressure.  The  jet  velocity  drops 
off  more  slowly  during  this  interval  than  during  the  buoyancy  dominated  motion. 

To  evaluate  the  jet  velocity  at  x9,  we  make  use  of  the  constancy  of  momentum  in 

* 11 

the  direction  normal  to  the  jet  axis,  that  is 


Pluml6lsin  “1  =P2um2  62  sin  "2 


(90) 


where  6 is  the  jet  width  and  a is  the  angle  between  the  jet  velocity  and  the  wind. 
Since  the  jet  is  thoroughly  mixed,  the  densities  Pl  and  p2  may  be  taken  as  equal. 
The  mixing  coefficient  c is  more  complex  here  than  in  the  case  of  a horizontal  jet 
since  the  angle  between  the  wind  and  jet  directions  is  no  longer  constant.  For 
simplicity,  we  shall  — following  Abramovich  — assume  that  the  value  of  c is  the 
same  as  for  a horizontal  jet.  The  value  of  *2  is  then  easily  calculated.  A numeri- 
cal check  of  Eq.  (90)  for  many  tests  show  that  u2  „ is  close  to  1/2  of  u2  , and. 
for  convenience,  we  shall  use  the  value  1/2. 

As  the  vertical  jet  interact  with  the  wind,  it  is  gradually  bent  over  toward  a 
horizontal  direction.  To  obtain  this  upper  trajectory,  we  employ  the  method 
developed  by  Shandorov  for  obtaining  the  path  of  a jet  in  a deflecting  flow.  12  On  the 


11.  Abramovich,  G.  N.  (1963)  The  Theory  of  Turbulent  Jets.  The  MIT  Press, 

Cambridge,  Mass.,  Chapter  12,  pp  547-548. 

12.  Abramovich,  G.N.  (1963)  The  Theory  of  Turbulent  Jets.  The  MIT  Press, 

Cambridge,  Mass.,  Chapter  12,  pp  541-556. 
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basis  of  many  experiments,  Shandorov  gives  an  empirical  equation  which,  for  the 
case  of  wind  perpendicular  to  the  jet,  may  be  written  as 


. . . .2.  55 

x2  ' x)  / ye  ~ y \ J_ 

26  \ 26  ) pe  u2 


(91) 


where  ug  is  the  initial  jet  velocity,  assumed  uniform  across  the  jet,  p the  jet 
density,  and  yg  the  initial  vertical  coordinate,  at  the  position  Xg.  Since  the  jet  is 
well  mixed  with  air,  we  may  take  pg  as  unity.  We  shall  assume  that  the  velocity 
profile  is  linear  at  Xg  and,  therefore,  calculate  the  effective  value  of  ug  by 


u 


2 

e 


where  u 


2 

m2 


is  taken  as  0.  5. 


(92) 


6.  DETERMINATION  OF  LIFT-OFF  POINT  AND  LENGTH  OF  JET  FROM 
EXPERIMENTAL  DATA 


In  correlating  the  lift-off  point  x^  against  the  experimental  data,  we  shall 

assume  that  x^  is  proportional  to  the  product  of  the  initial  jet  velocity  uq  and  the 

relative  velocity  ua<m  + 1)  and  inversely  proportional  to  the  temperature  excess 

ATm  and  write 
m 


x 

a 


u u (m  + 1) 
o a 

g y AT 
6 Jo  m 


(93) 


Using  Eq.  (3)  for  the  temperature  excess  ATm,  and  noting  that  m is  large  compared 
to  unity,  Eq,  (93)  becomes 

*a  ~ m2K  (94) 

where  K is  the  Froude  number  for  the  initial  flow.  Other  powers  of  the  initial  and 
relative  velocities  may  be  tried  in  Eq.  (93),  leading  to  different  combinations  of  m 
and  K in  Eq.  (94).  However,  the  relation  given  by  Eq.  (94)  appeared  to  yield  the 
best  results. 

The  results  of  the  correlation  are  shown  in  Figure  2 where  we  have  plotted  the 

1 1 2 

experimental  values  of  the  lift-off  point  x^  against  the  parameter  p = mK  ' „ 
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Figure  2.  Plot  of  Experimental  Values  of  Lift-Off 
Point  x Against  Parameter  p;  Solid  Curve  is  Power 
Law  Fit  to  Experimental  Points 


A fair  correlation  is  obtained,  but  not  as  good  as  that  obtained  for  the  coflowing 
jet.  The  data  has  been  fitted  with  a power  law  curve  with  the  result, 

x = 0.  66  p^‘  (95) 

a 

and  a correlation  coefficient  of  0.  79. 

Although  the  jet  length  x2  may  be  calculated  from  the  model,  as  shown  in  ref- 
erence 8,  it  is  desirable  to  be  able  to  obtain  x2  from  the  parameters  of  the  experi- 
ment. Since  the  jet  is  buoyant,  it  should  depend  upon  the  same  parameter  used  in 
determining  the  lift-off  point.  In  addition,  in  view  of  the  initial  heating  of  the  jet, 
x2  should  drop  off  inversely  to  the  inlet  jet  temperature  Tq.  However,  the  use  of 
Tq  in  the  correlation  did  not  result  in  any  visible  improvement  of  the  results.  We 
have,  therefore  — for  convenience  — utilized  the  same  parameter  p for  jet  length  as 
previously  used  for  lift-off  point.  The  results  of  the  correlation  are  shown  in  Fig- 
ure 3 where  x2  is  plotted  against  p.  The  amount  of  scatter  is  about  the  same  as 
that  obtained  for  x^.  A power  law  curve  has  been  fitted  to  the  data  yielding 
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0 71 

x2  = 4.02  pu‘  (96)  4 

| 

with  a correlation  coefficient  of  0.  83. 


Figure  3.  Plot  of  Experimental  Values 
of  Jet  Length  X2  Against  Parameter  p; 

Solid  Curve  is  Power  Law  Fit  to  Exper-  < 

imental  Points 


7.  RESULTS  AND  DISCUSSION 


The  calculated  trajectories  are  shown  in  Figures  4-9  along  with  the  correspond- 
ing experimental  curves.  These  results  cover  most  of  the  range  of  windspeeds, 
initial  jet  velocities,  and  initial  jet  temperatures  encountered  in  the  tests.  The 
calculated  results  are  in  fair  agreement  with  the  experimental  curves,  but  the 
lower  portion  of  the  calculated  trajectory  generally  tends  to  be  somewhat  higher 
tnan  the  corresponding  experimental  result.  The  lower  buoyancy  observed  experi- 
mentally may  be  due  to  the  fact  that  when  the  jets  merge,  the  assumed  equivalent 
planar  jet  contains  regions  of  cooler  air  between  the  jets  which  help  lower  the  total 
buoyant  force  of  the  jet.  In  addition,  when  the  jet  detaches  from  the  ground  some 
jet  air  may  become  entrained  with  the  ambient  air  going  under  the  jet,  resulting  in 


A. 
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Figure  6.  Comparison  of  Calculated  Trajectories,  for  Several  Tests, 
with  Corresponding  Experimental  Curves,  Test  5-4 


Figure  7.  Comparison  of  Calculated  Trajectories,  for  Several  Tests, 
with  Corresponding  Experimental  Curves,  Test  5-19 
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Figure  8.  Comparison  of  Calculated  Trajectories,  for  Several  Tests, 
with  Corresponding  Experimental  Curves,  Test  6-2 


Figure  9.  Comparison  of  Calculated  Trajectories,  for  Several  Tests, 
with  Corresponding  Experimental  Curves,  Test  8-10A 
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lower  average  buoyant  force.  An  indication  of  these  effects  is  shown  in  Test 
5-19  (Figure  7)  for  which  the  calculated  curve  shows  a steady  moderate  rise  after 
lift-off  in  contradistinction  to  the  experimental  trajectory  which  remains  low  a 
considerable  distance  before  starting  to  rise  rapidly. 

In  several  of  the  tests,  5-4  and  6-2  (Figures  6 and  8),  the  calculated  curves 
appear  to  rise  at  a proper  rate  but,  because  of  an  incorrect  jet  length  or  lift-off 
point,  deviate  moderately  from  the  experimental  trajectory.  This  behavior  is  due 
to  the  scatter  of  the  lift-off  points  and  jet  length  points  (Figures  2 and  3),  resulting, 
in  several  cases,  in  appreciable  deviation  from  the  fitted  curves. 

In  view  of  the  simplified  model  utilized  for  the  upper  portion  of  the  trajectory, 
the  agreement  between  the  calculated  and  experimental  results  is  quite  reasonable. 
Additional  refinement  may  be  incorporated  into  the  model  described  in  Section  5, 
that  is,  a more  detailed  calculation  of  the  turbulent  mixing  and  bending  of  the  curve 
as  it  approaches  the  vertical.  It  is  felt,  however,  that  such  refinement  should  be 
incorporated  into  a more  detailed  description  and  development  of  the  entire  model. 


8.  SUMMARY  AND  CONCLUSIONS 

A model,  previously  developed  for  obtaining  the  dynamic  properties  of  a heated 
turbulent  counterflowing  jet  in  the  absence  of  buoyancy,  has  been  utilized  to  obtain 
the  lower  portion  of  the  jet  trajectory  when  the  buoyancy  is  taken  into  account. 

The  upper  portion  of  the  trajectory  was  obtained  from  a model  which  yields  the  de- 
flection of  a jet  in  a cross  wind.  In  general,  the  calculated  trajectories  are  in  fair 
to  good  agreement  with  the  corresponding  experimental  curves,  the  calculated  re- 
sults tending  to  be  somewhat  higher  than  those  obtained  experimentally. 

A scaling  law  has  been  obtained  for  the  counterflowing  jet,  which  indicates  that 
the  scaling  depends  principally  upon  the  Froude  number  referred  to  initial  jet 
velocity  and  excess  temperature  and  very  weakly  upon  initial  jet  temperature.  This 
result  is  essentially  the  same  as  that  obtained  for  the  coflowing  jet  for  the  case  of 
small  values  of  wind  speed  relative  to  initial  jet  velocity. 
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List  of  Symbols 


b distance  from  outside  wall  to  jet  axis 

c jet  thickness  coefficient 

H end  of  initial  section 

K Froude  number 

m wind  speed  parameter  = UD/Ua 

p jet  parameter  = 

T jet  temperature 

T&  ambient  temperature 

T temperature  on  axis 

m 

AT  temperature  excess  over  ambient  (T  - T&) 

u jet  velocity 

um  velocity  on  axis 

V vertical  velocity 

x horizontal  position  along  jet 

x^  end  of  initial  section 

x^  end  of  Zone  I 

Xj  end  of  Zone  II 
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lift-off  position 
vertical  position 
jet  half-width 

zero  velocity  surface  in  main  section 

surface  separating  perturbed  and  unperturbed  flows 

angle  between  jet  and  wind 

gas  density 

jet  width 


Subscripts 

a 

ambient 

c 

region  above  y^  surface 

m 

on  jet  axis 

o 

initial  jet  position 

a 

lift-off  position 

e 

initial  position 

y 

y0 

*2 

*3 

a 

P 
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